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Abstract: We discuss the lower Higgs boson mass bounds which come from the absolute 
stability of the Standard Model (SM) vacuum and from the Higgs inflation, as well as the 
prediction of the Higgs boson mass coming from asymptotic safety of the SM. We account 
for the 3-loop renormalization group evolution of the couplings of the Standard Model and 
for a part of two-loop corrections that involve the QCD coupling a s to initial conditions 
for their running. This is one step above the current state of the art procedure ( "one-loop 
matching-two-loop running"). This results in reduction of the theoretical uncertainties 
in the Higgs boson mass bounds and predictions, associated with the Standard Model 
physics, to 1-2 GeV. We find that with the account of existing experimental uncertainties 
in the mass of the top quark and a s (taken at 2a level) the bound reads Mu > M m ; n 
(equality corresponds to the asymptotic safety prediction), where M m \ n = 129 ± 6 GeV. 
We argue that the discovery of the SM Higgs boson in this range would be in agreement 
with the hypothesis of the absence of new energy scales between the Fermi and Planck 
scales, whereas the coincidence of Mh with M m ; n would suggest that the electroweak scale 
is determined by Planck physics. In order to clarify the relation between the Fermi and 
Planck scale a construction of an electron-positron or muon collider with a center of mass 
energy ~ 200 + 200 GeV (Higgs and t-quark factory) would be needed. 
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1 Introduction 

The mass Mh of the Higgs boson in the Standard Model is an important indicator of the 
presence of new energy scales in particle physics. It is well known that if M™^ a < Mh < 
M max dau then the SM is a consistent effective field theory all the way from the Fermi scale 
up to the (reduced) Planck scale Mp = 2.44 x 10 18 GeV. The upper limit comes from the 
requirement that the Landau pole in the scalar self-coupling 1 must not appear at energies 
below Mp [1-3]. The lower limit comes from the requirement of the stability of the SM 
vacuum against tunneling to the states with the Higgs field <\> exceeding substantially the 
electroweak value 250 GeV [4-6] (see Fig. 1). 
The estimates of M^ dau 

give a number around 175 GeV [1-3, 7] which is in the Mh 
range excluded (at least in the range 129 — 525 GeV) by the searches for the SM Higgs 

lr To be more precise, the scalar self-coupling is infinite in the one-loop approximation only. If higher 
order terms are included, it may not become infinite, but move away from the region of the weak coupling. 
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boson at the LHC and Tevatron [8, 9]. In other words, we already know that the SM is a 
weakly coupled theory up to the Planck scale. 

One can distinguish between two types of the stability bounds. If Mu > M^ blllty , 
the electroweak vacuum is absolutely stable, whereas if M™f n ta < M H < M^? bflity , then 
it is metastable with the life-time exceeding that of the Universe. Numerically, M™^ a ~ 
111 GeV [10]. The existence of the Higgs boson with the mass smaller than M™-? a would 
provide an undisputed argument in favor of existence of new physics between the Fermi 
and Planck scale. However, already since LEP we know that this is not the case. 

The Higgs mass M^ blllty is not at all special from the point of view of the validity 
of the SM up to the Planck scale. The value of M^ blllty , however, plays a crucial role 
if the Standard Model is embedded to a bigger picture which includes gravity. First, 
only if Mn > Mj^ >lllty , the Higgs boson of the SM can play the role of the inflaton and 
make the Universe flat, homogeneous and isotropic, and produce a necessary spectrum 
of perturbations needed for structure formation [11, 12]. Second, Mjj = M^ blllty is a 
prediction of asymptotically safe scenario for the SM [13], making it consistent up to 
arbitrary large scale. 

Thus, we will focus on the upgrade of existing computations of M^ blllty and on the 
discussion of the significance of the relation between the Higgs boson (to be discovered yet) 
mass M H and M^ bmty for beyond the SM (BSM) physics. 

The computation of M^' i '^ blllty has been already done in a large number of papers 
[10, 12, 14-17]. It is divided into two parts. The first one is the determination of the 
MS parameters from the physical observables and the second one is the renormalization 
group running of the MS constants from the electroweak to a high energy scale. The 
most advanced recent works [10, 12, 17, 18] use the so-called "one-loop-matching-two-loop- 
running" procedure. It can determine the Higgs boson mass bounds with the theoretical 
accuracy of 2 — 5 GeV (see the discussion of uncertainties in [12] and below). Meanwhile, 
the most important terms in the 3-loop running of the gauge and Higgs coupling constants 
were computed in [19, 20] (we thank K. Chetyrkin and M. Zoller for sharing these results 
with us prior to publication). The present work accounts for 0{aa s ) corrections in the 
MS-pole matching procedure, which were not known previously. This allows us to decrease 
the theoretical uncertainties in the Higgs boson mass prediction/bounds, associated with 
the SM physics down to 1-2 GeV. This is a new result, based on a superior partial "two- 
loop-matching-three-loop-running" procedure. These findings are described in Section 2.2. 
We will see that the experimental errors in the mass of the top-quark and in the value 
of the strong coupling constant are too large to settle up the question of the stability of 
the electroweak vacuum, even if the LHC will confirm the evidence for the Higgs signal 
presented by the ATLAS and CMS collaborations [8, 9] in the region M H = 124-126 GeV. 

In Section 3 we will discuss the significance of the relationship between the true Higgs 
boson mass Mjj and M^ blllty for BSM physics. We will argue that if Mjj = M^ blllty then 
the electroweak symmetry breaking is likely to be determined by Planck physics and that 
this would indicate an absence of new energy scales between the Fermi and gravitational 
scales. We will also address here the significance of M^? lity for the SM with gravity 
included. Of course, this can only be done under certain assumptions. Specifically, we 
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Figure 1. Higgs self-coupling in the SM as a function of the energy scale. The top plot depicts 
possible behaviors for the whole Higgs boson mass range — Landau pole, stable, or unstable elec- 
troweak vacuum. The lower plots show detailed behavior for low Higgs boson masses, with dashed 
(dotted) line corresponding to the experimental uncertainty in the top mass M t (strong coupling 
constant a s ), and the shaded yellow (pink) regions correspond to the total experimental error and 
theoretical uncertainty, with the latter estimated as 1.2 GeV (2.5 GeV), see section 2 for detailed 
discussion. 



will discuss the non-minimal coupling of the Higgs field to the Ricci scalar (relevant for 
Higgs-inflation [11, 12, 21]) and the asymptotic safety scenario for the SM [13]. 

In Section 4 we present our conclusions. We will argue that if only the Higgs boson with 
the mass around M^ 1 blllty and nothing else will be found at the LHC, the next step in high 
energy physics should be the construction a new electron-positron (or muon) collider — the 
Higgs and t-factory. It will not only be able to investigate in detail the Higgs and top 
physics, but also elucidate the possible connection of the Fermi and Planck scales. 

Appendix A contains a full set of formulas required for the determination of the MS 
coupling constants from the pole masses of the SM particles, including the corrections of 
the orders of up to O(af), 0(a), and 0(aa s ). The computer code for the matching is 
made publicly available at http : //www. inr . ac . ru/~f edor/SM/. 
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Fermi Planck Fermi Planck 

Figure 2. Schematic depiction of the SM effective potential V for the Higgs field for Mh > 
Mobility ( i eft) and Mh < M ™ y (right). 

2 The stability bound 

The stability bound will be found in the "canonical" SM, without any new degrees of 
freedom or any extra higher dimensional operators added, see Fig. 2. 

2.1 The benchmark mass 

It will be convenient for computations to introduce yet another parameter, "benchmark 
mass", which we will call M m [ n (without any superscript). Suppose that all parameters of 
the SM, except for the Higgs boson mass, are exactly known. Then M m i n , together with 
the normalisation point /j,q, are found from the solution of two equations: 

A(/x ) = 0, (3 x (\(po)) = 0, (2.1) 

where j3\ is the /3-function governing the renormalisation group (RG) running of A. Here 
we define all the couplings of the SM in the MS renormalisation scheme which is used 
de-facto in the most of the higher-loop computations. Clearly, if any other renormalization 
scheme is used, the equations A = j3\ = will give another benchmark mass, since the 
definition of all the couplings are scheme dependent. 

The procedure of computing M m ; n is very clean and transparent. Take the standard 
MS definition of all coupling constants of the SM, fix all of them at the Fermi scale given 
the experimentally known parameters such as the mass of the top quark, QCD coupling, 
etc., and consider the running Higgs self-coupling X(/j,) depending on the standard t'Hooft- 
Veltman parameter fj,. Then, adjust M m \ n in such a way that equations (2.1) are satisfied 
at some /j,q. 

For the stability bound one should find the effective potential V((f>) and solve the 
equations 

V(cj> SM ) = V{ct> 1 ), V'^sm) = V'fai) = 0, (2.2) 

where ())sm corresponds to the SM Higgs vacuum, and (pi correspond to the extra vacuum 
states at large values of the scalar field. Though the effective potential and the field <j> 
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are both gauge and scheme dependent, the solution for the Higgs boson mass to these 
equations is gauge and scheme invariant. 

In fact, M^j* blllty is very close to M m - m . Numerically, the difference between them 
is much smaller, than the current theoretical and experimental precisions for M m i n , see 
below. The following well known argument explains why this is the case. The RG improved 
effective potential for large 4> can be written as [15, 16, 22] 



V(</>) oc A(</>)</> 4 [l + O MMi/Mj))] , (2.3) 



where a is here the common name for all the SM coupling constants (which are rather 
small at Planck scale), and Mi are the masses of different particles in the background of 
the Higgs field. If 0(a) corrections are neglected, the equations (2.2) coincide with (2.1), 
meaning that M m ; n ~ M^^ lhty . The numerical evaluation for one loop effective potential 
gives Am stablllty = M^^ lllty — M m i n ~ — 0.15GeV, which can be neglected in view of 
uncertainties discussed below. 

Note that in many papers the stability bound is shown as a function of the cutoff scale 
A (the energy scale up to which the SM can be considered as a valid effective field theory). 
It is required that V(<p) > V(4>sm) for all cj) < A. This can be reformulated as X(fi) > 
for all (i < A with pretty good accuracy. Interestingly, if A = Mp, this bound is very 
close to the stability bound following from eq. (2.2), having nothing to do with the Planck 
scale (see also below). Note also that the uncertainties in experimental determinations of 
M t and a s together with theoretical uncertainties, described in the next section, lead to 
significant changes in the scale A. Fig. 1 illustrates that for Higgs boson masses 124 — 127 
GeV this scale may vary from 10 s GeV up to infinity within currently available precisions. 

2.2 Value of M min 

The state of art computation of M m \ a contained up to now the so called one-loop MS-pole 
matching, relating the experimentally measured physical parameters to the parameters of 
the SM in the MS subtraction scheme (to be more precise, the two-loop a s corrections to 
the top pole mass-MS mass relation has been included [10]). Then the results of the first 
step are plugged into two-loop RG equations and solved numerically. 

Before discussing the upgrade of the one-loop-matching-two-loop-running procedure, 
we will remind of the results already known and their uncertainties. We will make use of 
our computations of M m \ n presented in [12]. 2 A somewhat later papers [17, 18, 23] contains 
exactly the same numbers for M^^ lhty (note, however, that the theoretical uncertainties 
were not discussed in [23]). See also earlier computations in [7, 10, 14-16, 18]. 

In [12] we found: 



M„ 



M t - 171.2 GeV , a s - 0.1176 

126.3 + — x 4.1 x 1.5 

2.1 GeV 0.002 



GeV, (2.4) 



2 The main interest in this paper was the lower bound on the Higgs boson mass in the Higgs-inflation, 
see below. However, M m i n has been estimated as well as a by-product of the computation. 
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Source of uncertainty 


Nature of estimate 


A t v lonr M rr ,jn, GeV 


*i-ioop niauciiiiig A 


"^£}T1 CI T IT n TT T T/~» II 

oeiioiLi viiy lo fj, 


1 7 


2-loop matching y t 


Sensitivity to /i 


0.6 


3-loop a s to yt 


known 


1.4 


4-loop a s to yt 


educated guess [27, 28] 


0.4 


confinement, j/t 


educated guess ~ ^-qcd 


n k 
u.o 


3-loop running M\\r — y Mp 


educated guess 


0.8 


total uncertainty 


sum of squares 


2.5 


total uncertainty 


linear sum 


5.4 


Corrections to [12] 




AM min , GeV 


Typos in the code used in [12] 


error 


+0.2 


Extra 5^ m (A. 5) of [10] 


error 


+0.4 


"Exact" formula instead of 






approximation (2.20) in [32] 


clarification 


+0.1 


Total correction to (7.1) of [12] 




+0.7 


Total shift to be applied to (7.1) 


of [12] for comparison 


+0.7 



Table 1. Theoretical uncertainties and mistakes in the M m i n evaluation in [12]. 

and estimated the theoretical uncertainties as summarized in Table 1 (see also [17] ). 3 While 
repeating this analysis we found some numerical errors which are given at the bottom 
section of this table (see a detailed discussion below). In total, they shift the value given 
in eq. (2.4) up by 0.7 GeV. As for uncertainties, they were estimated as follows. The 
one- loop matching formulas can be used directly at /i = rrit, or at some other energy 
scale, e.g. at /i = Mz, and then the coupling constants at m t can be derived with the 
use of RG running. The difference in procedures gives an estimate of two-loop effects in 
the matching procedure. This is presented by the first two lines in Table 1 (in fact, we 
underestimated before the uncertainty from A matching — previously we had here 1.2 GeV 
and now 1.7 GeV). The next two lines are associated with 3 and 4-loop corrections to 
the top Yukawa yt. The 3-loop corrections were computed in [24-26] and the four- loop 
a s contribution to the top mass was guessed to be of the order 5yt(mt) / yt{fnt) — 0.001 in 
[27, 28]. The non-perturbative QCD effects in the top pole mass-MS mass matching are 
expected to be at the same level [29-31]. For 3-loop running we put the typical coefficients 
in front of the largest couplings a s and yt- If these uncertainties are not correlated and 
can be summed up in squares, the theoretical uncertainty is 2.5 GeV. If they are summed 
up linearly, then the theoretical error can be as large as ~ 5 GeV. 

Now, this computation can be considerably improved. First, in [19] the 3-loop cor- 

3 The reader not interested in details of the comparison with the previous results can skip this discussion 
and restart at the following paragraph. 
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Contribution AM m i n , GeV 

Three loop beta functions -0.23 
5y t oc 0{a 3 s ) -1.15 
5yt oc 0(aa s ) -0.13 
5X oc 0(aa s ) 0.62 



Table 2. Contributions to the value of the M m i n . 

rections to the running of all gauge couplings has been calculated. Second, in [20] the 
leading contributions (containing the top Yukawa and a s ) to the running of the top quark 
Yukawa and the Higgs boson self coupling have been determined. This removes the uncer- 
tainty related to 3-loop RG running. In addition, in the present paper, we determine the 
two- loop corrections of the order of 0{aa s ) to the matching of the pole masses and the 
top quark Yukawa and Higgs boson self coupling constants. Also, the known [24-26] three 
loop QCD correction to the top quark mass relation of the order O(af) can be included 
(previously it had been used for estimates of uncertainties) . All this considerably decreases 
the theoretical uncertainties in M m i n . 

The individual contributions of the various new corrections on top of the previous 
result are summarized in the Table 2. It is clearly seen that there are two new significant 
contributions — one is the three-loop pure QCD correction to the top quark mass [24-26], 
and another is the two loop correction 0(aa s ) to the Higgs boson mass, found in the 
present paper. Together the new contributions sum to the overall shift of the previous 
prediction [12] by —0.89 GeV, giving the result 4 



M„ 



M t - 172.9 GeV a s - 0.1184 



GeV. (2.5) 



The new result (2.5) is less than 0.2 GeV away from the old one (2.4) if the same central 
values for Mt and a s are inserted. This coincidence is the result of some magic. In the old 
evaluation several mistakes were present, summarized in Table 1. The largest one was the 
double counting of <5^ ED in (A. 5) of [10], as compared to the original result [32]. Also, there 
were minor typos in the computer code for the matching of the Higgs coupling constant, 
and finally there was a small correction coming from the use of an approximate rather 
than exact one loop formula for O(a) corrections from [32]. These corrections add 0.7 GeV 
to the original number in [12]. By chance this almost exactly canceled the —0.89 GeV 
contribution from the higher loops, Table 2, nearly leading to a coincidence of (2.5) and 
(2.4). 

Table 3 summarizes the uncertainties in the new computation. It contains fewer lines. 
Now we can ignore safely the error from higher order (4-loop) RG corrections for the 
running up to the Planck scale. The first two lines were derived in the same manner as 



4 Note, that this value is the benchmark mass defined by eq. (2.1). The mass, corresponding to the 
metastability bound (2.2) is ~ 0.1-0.2 GeV smaller. 
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Source of uncertainty 



Nature of estimate A t h e0 r-^min, GeV 



3-loop matching A 

3- loop matching y t 

4- loop a s to yt 
confinement, yt 

4-loop running Myy Mp 



Sensitivity to fj, 
Sensitivity to \x 
educated guess [27, 28] 
educated guess ~ Aqcd 
educated guess 



1.0 
0.2 
0.4 
0.5 
< 0.2 



total uncertainty 
total uncertainty 



sum of squares 
linear sum 



2.3 



1.2 



Table 3. Theoretical uncertainties in the present M min evaluation. 



previously. For the Higgs boson self-coupling we can use the matching formulas (A. 35) 
to get the value of X(/x) at scale \i = M t directly, or to get the value \{Mz) and then 
evolve the constants to the scale fx = Mt with the RG equations. The obtained difference 
5X(Mt) I 'X(Mt) ~ 0.016 corresponds to the error 5m ~ 1.0 GeV. A similar procedure of 
comparing evolution between Mt and Mz using RG equations and direct matching formulas 
to the order 0(al,a,aa s ) leads for the chang in the top quark Yukawa bytjyt ~ 0.0005, 
leading to 5M ~ 0.2 GeV. Note, however, that strictly speaking this test verifies the error of 
the fi dependent terms in the matching formulas, while the constant ones may lead to larger 
contributions. We also do not estimate now the contributions of the order 0(a 2 ), where 
formal order in a may correspond to yf . Thus, this estimate should be better considered 
as a lower estimate of the error. The 4-loop matching and confinement contributions are 
the same as before. 

As an indication of the dependence on the matching point we present Fig. 3, where the 
reference Higgs boson mass M m ; n was obtained using the matching formulas at scale /io 
varying between the Z-boson and top quark masses. One can see that the overall change 
of the Higgs boson mass is about GeV. 

If we assume that these uncertainties are not correlated and symmetric we get a the- 
oretical error in the determination of the critical Higgs boson mass 5m^ e0 r — 1-2 GeV. If 
they are summed up linearly, we get an error of 2.4 GeV. We leave it to the reader to de- 
cide which estimate of the uncertainties is more adequate. The precision of the theoretical 
value of M m i n can be further increased by computing the 0(a 2 ) two- loop corrections to the 
matching procedure. Numerically, the most important terms are those when a corresponds 
to y 2 and A. 

The result (2.5) is visualized by Fig. 4. The experimentally allowed regions for the 
top mass Mt and strong coupling a s are adopted PDG 2010 edition [33]. 5 On top of these 
allowed regions the bands corresponding to the reference values of the Higgs boson mass 
M m i n being equal to 124, 125, 126, 127 GeV are shown, with the dashed and dotted lines 

J Note however, that the current experimental error estimate is based on averaging over different exper- 
imental approaches. In some methods quite a different central values are obtained. See e.g. [34-36] about 
a s determination and [37-39] about M t . 
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Figure 3. The dependence of the reference Higgs boson mass M m i n on the matching scale /xo (the 
MS constants are obtained by matching formulas at scale /xn and then used for the solution of the 
equations (2.1)). The thick solid line corresponds to the full matching formulas A ~ 0(a,aa s ), 
yt ~ 0(or s , a, aa s ); the thin lines correspond to using matching formulas of lower order. The thick 
dashed line corresponds to using additionally the two loop electroweak contributions to the higgs 
coupling constant in the gauge-less limit, eq. (48) of [112], see discussion in "Note added". Here 
M t = 172.9 GeV and a s = 0.1184. 

corresponding to quadratically or linearly added estimates of theoretical uncertainties. 

One can see that the accuracy of theoretical computations and of the experimental 
measurements of the top and the Higgs boson masses does not allow yet to conclude with 
confidence whether the discovery of the Higgs boson with the mass 124 — 127 GeV would 
indicate stability or metastability of the SM vacuum. All these reference values of Higgs 
masses are compatible within 2a with current observations. 

3 M min and BSM physics 

Our definition of the "benchmark" Higgs boson mass consists of the solution of the two 
equations (2.1) and gives, in addition to M m i n , the value of the scale no at which the 
scalar self-coupling and its /3-function vanish simultaneously. The central value for is 
2.9 x 10 18 GeV and is quite stable if mt and a s are varied in their confidence intervals 
(see Fig. 5). One can see that there is a remarkable coincidence between an d the 
(reduced) Planck scale Mp = 2.44 x 10 18 GeV. The physics input in the computation 
of (io includes the parameters of the SM only, while the result gives the gravity scale. 
A possible explanation may be related to the asymptotic safety of the SM, see [13] and 
below. 6 It remains to be seen if this is just the random play of the numbers or a profound 
indication that the electroweak symmetry breaking is related to Planck physics. If real, 
this coincidence indicates that there should be no new energy scales between the Planck 
and Fermi scales, as they would remove this coincidence unless some conspiracy is taking 
place. 

Of course, this is not the only way to predict the Higgs boson mass close to M m i n . Cf. the "multiple 
point principle" of [40, 41], requiring the degeneracy between the SM vacuum and an extra one appearing 
at the Planck scale, the inflation from false vacuum decay [42-45], etc. 
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Figure 4. The values of the strong coupling constant a s and top mass M t corresponding to several 
minimal Higgs boson mass M m i n . The 68% and 95% experimentally allowed regions for a s and M t 
are given by shaded areas. The dashed (dotted) lines correspond to 1.2 GeV (2.45 GeV) uncertainty 
in the M m - ln theoretical determination. 



We will discuss below two possible minimal embeddings of the SM to the theory of 
gravity and discuss the significance of M m ; n in them. 

3.1 Asymptotic safety 

The asymptotic safety of the SM [13], associated with the asymptotic safety of gravity [46], 
is strongly related to the value of the Higgs boson mass. Though General Relativity is non- 
renormalizable by perturbative methods, it may exist as a field theory non-perturbatively, 
exhibiting a non-trivial ultraviolet fixed point (for a review see [47]). If true, all other 
coupling of the SM (including the Higgs self-interaction) should exhibit an asymptotically 
safe behaviors with the gravity contribution to the renormalisation group running included. 

The prediction of the Higgs boson mass from the requirement of asymptotic safety 
of the SM is found as follows [13]. Consider the SM running of the coupling constants 
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Figure 5. The scale fio (solution of (2.1)) depending on the top mass M t . The dashed lines 
correspond to la uncertainty in the a s . The yellow shaded region corresponds to adding the a s 
experimental error and the theoretical uncertainty in the matching of the top Yukawa yt and top 
pole mass. 



and add to the /3-functions extra terms coming from gravity, deriving their structure from 
dimensional analysis: 

2 

Ph ~~ 87r M|(/i) [6A) 

where a±, a,2, 03, a y , and a\ are some constants (anomalous dimensions) corresponding 
to the gauge couplings of the SM g, g', g s , the top Yukawa coupling y t , and the Higgs 
self-coupling A. In addition, 

M|( M ) ~ Mp + (3.2) 

is the running Planck mass with £0 ~ 0.024 following from numerical solutions of functional 
RG equations [48-50]. Now, require that the solution for all coupling constants is finite for 
all fi and that A is always positive. The SM can only be asymptotically safe if a\, 02, 03, a y 
are all negative, leading to asymptotically safe behavior of the gauge and Yukawa couplings. 
For a\ < we are getting the interval of admissible Higgs boson masses, M^^ ty < Mjj < 
Mrm« ty - However, if a\ > 0, as follows from computations of [49, 50], only one value of 
the Higgs boson mass Mh = M^^ ty leads to asymptotically safe behavior of A. As is 
explained in [13], this behavior is only possible provided X(Mp) and f3\(X(Mp)) 0. 
And, due to miraculous coincidence of /io and Mp, the difference Am safety = M^^ ty — M m ; n 
is extremely small, of the order 0.1 GeV. The evolution of the Higgs self-coupling for the 
case of ah < is shown in Fig. 6, and for the case ah > in Fig. 7. 

In fact, in the discussion of the asymptotic safety of the SM one can consider a more 
general situation, replacing the Planck mass in eq. (3.2) by some cutoff scale A = nMp. 
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Figure 6. Schematic depiction of the behavior of the scalar self-coupling if ah < for M^^ ty < 
Mh < Ai^fax ty (left) and Mjj < M^^ ty (right). In both cases gravity leads to asymptotically free 
behavior of the scalar self-coupling. Negative A lead to instability and thus excluded. 
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Figure 7. Schematic depiction of the behavior of the scalar self-coupling if ah > for Mh > 
M^;^ ty , leading to Landau-pole behavior (left), Mh > M^?® ty , leading to instability (right) and 
Mh = A^n ty > asymptotically safe behavior (middle). Only this choice is admissible. 



Indeed, if the Higgs field has non-minimal coupling with gravity (see below), the behavior 
of the SM coupling may start to change at energies smaller than Mp by a factor l/£, 
leading to an expectation for the range of k as 1/^ < /t < 1. Still, the difference between 
M m ; n and M^^ ty remains small even for k ~ 10~ 4 , where M^° ty ~ 128.4 GeV, making the 
prediction Mh — M m \ a sufficiently stable against specific details of Planck physics within 
the asymptotic safety scenario. 

3.2 Mmin and cosmology 

It is important to note that if the mass of the Higgs boson is smaller than the stability 
bound M m i n , this does not invalidate the SM. Indeed, if the life-time of the metastable 
SM vacuum exceeds the age of the Universe (this is the case when Mh > M me t a , with 
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M meta — 111 GeV [10]) then finding a Higgs boson in the mass interval M meta < Mjj < 
M m ; n would simply mean that we live in the metastable state with a very long lifetime. 
Of course, if the Higgs boson were discovered with a mass below M meta , this would prove 
that there must be new physics between the Fermi and Planck scales, stabilizing the SM 
vacuum state. However, the latest LEP results, confirmed recently by LHC, tell us that in 
fact Mh > -^metai and, therefore, that the presence of a new energy scale is not required, 
if only the metastability argument is used. 

The bound Mh > -Mmeta can be strengthened if thermal cosmological evolution is con- 
sidered [10]. After inflation the universe should find itself in the vicinity of the SM vacuum 
and stay there till present. As the probability of the vacuum decay is temperature depen- 
dent, the improved Higgs boson mass bound is controlled by the reheating temperature 
after inflation (or maximal temperature of the Big Bang). The latter is model dependent, 
leading to the impossibility to get a robust bound much better than M meta . For example, 
in R? inflation [51, 52] the reheating temperature is rather low, T ~ 10 9 GeV [52], leading 
to the lower bound 116 GeV [53] on the Higgs boson mass, which exceeds M meta only by 
4 GeV. 

However, if no new degrees of freedom besides those already present in the SM are 
introduced and the Higgs boson plays the role of inflaton, the bound Mjj > M m j n reappears, 
as is discussed below. 

3.3 Higgs inflation 

The inclusion of a non-minimal interaction of the Higgs field with gravity, given by the 
Lagrangian ^|^| 2 i?, where R is the Ricci scalar, changes drastically the behavior of the 
Higgs potential in the region of large Higgs fields 4> > ^inflation — ^p/VC [H]- Basically, 
the potential becomes flat at 4> > Mi n fl ation , keeping the value it acquired at 4> ~ Mp/y/£. 
This feature leads to a possibility of Higgs-inflation: if the parameter £ is sufficiently large, 
700 < £ < 10 5 , [12] the Higgs boson of the SM can make the Universe flat, homogeneous 
and isotropic, and can produce the necessary spectrum of primordial fluctuations. The 
possibility of the Higgs inflation is also strongly related to the value of the Higgs boson 
mass: the successful inflation can only occur if M™ m atlon < M H < M™ a atlon . The upper 
limit M™ a atlon comes from the requirement of the validity of the SM up to the inflation 
scale M m fl at i on . Near M™ m atlon the behavior of the effective potential in the Einstein frame 
changes as shown in Fig. 8: if Mg < M™ m atlon the "bump" in the Higgs potential prevents 
the system to go to the SM vacuum state. As in the previous case, these bounds can be 
formulated with the use of the Higgs self-coupling A. Basically, it must be perturbative 
and positive for all energy scales below M in fl ation . Though any Higgs boson mass in the 
interval M™ m ation < M H < M™ a ation can lead to successful inflation, the value M^ m ation 
is somewhat special. For the lower part of the admitted interval the value of the non- 
minimal coupling £ reaches its minimal value £ ~ 700, extending the region of applicability 
of perturbation theory [12, 54, 55]. 

The computation of the lower bound on the Higgs boson mass from inflation is more 
complicated. It is described in detail in [12, 21]. Basically, one has to compute the Higgs 
potential in the chiral electroweak theory associated with large values of the Higgs field 
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Figure 8. Schematic depiction of the effective potential V for the Higgs field in the Higgs- 



inflationary theory in the Einstein frame for Mjj > m 



inflation 



(left) and Mh < m 



inflation 



(right). 



and find when the slow-roll inflation in this potential can give the large-scale perturbations 
observed by the COBE satellite. The outcome of these computations, however, can be 
formulated in quite simple terms: for inflationary bound find M^^ tlon from the condition 
A(/i) > for all \x < M m fl ation [12]. A priori, the inflationary bound could have been very 
different from M m ; n and thus from M^? 1 ty . Indeed, both M m \ n and M ^ blllty know nothing 
about the Planck scale and are defined entirely within the SM, whereas the inflationary 
bound does use Mp. However, the remarkable numerical coincidence, between no and Mp, 
makes M m [ n and M™? atlon practically the same. The coupling constant A evolves very 
slowly near the Planck scale, so that the regions for the Higgs boson mass following from 
the conditions A(/z) > for \x < Mp and fi < M m fl at i on are almost identical. This leads to 
the result that Am inflation = M^ ation - M min ~ -0.1 - 0.2 GeV. This number is derived 
within the SM without addition of any higher dimensional operators. 

One must note, that simple scale analysis leads to the unitarity violation in the Higgs 
inflation below the Planck energy scale [56-59]. This means, that calculations in the 
Higgs inflationary models should be done with some additional assumptions about the high 
energy physics, formulated in [54], specifically the approximate scale invariance at high field 
backgrounds. As is explained in [54], adding to the SM higher-dimensional operators with 
a Higgs-field dependent cutoff modifies the lower bound on the Higgs boson mass in Higgs 
inflation. If these operators are coming with "natural" power counting coefficients (for 
exact definition see [54]) the sensitivity of the Higgs boson mass bound to unknown details 
of ultraviolet physics is rather small AM* 1 ™ ~ 0.6 GeV [54]. At the same time, it is 
certainly not excluded that the change of M^^ tlon can be larger. 

4 Conclusions 



If the SM Higgs boson will be discovered at LHC in the remaining mass interval 115.5 < 
Mh < 127 GeV not excluded at 95% [8, 9], there is no necessity for a new energy scale 
between the Fermi and Planck scales. The EW theory remains in a weakly coupled region 
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all the way up to Mp, whereas the SM vacuum state lives longer than the age of the 
Universe. If the SM Higgs boson mass will be found to coincide with M m i n given by (2.5), 
this would put a strong argument in favor of the absence of such a scale and indicate that 
the electroweak symmetry breaking may be associated with the physics at the Planck scale. 

The experimental precision in the Higgs boson mass measurements at the LHC can 
eventually reach 200 MeV and thus be much smaller than the present theoretical (~ 
1-2 GeV) and experimental (~ 5GeV, 2a) uncertainties in determination of M m ; n . The 
largest uncertainty comes from the measurement of the mass of the top quark. It does 
not look likely that the LHC will substantially reduce the error in the top quark mass 
determination. Therefore, to clarify the relation between the Fermi and Planck scales 
a construction of an electron-positron or muon collider with a center-of-mass energy of 
~ 200 + 200 GeV (Higgs and t-quark factory) would be needed. This would be decisive 
for setting up the question about the necessity for a new energy scale besides the two ones 
already known — the Fermi and the Planck scales. In addition, this will allow to study in 
detail the properties of the two heaviest particles of the Standard Model, potentially most 
sensitive to any types on new physics. 

Surely, even if the SM is a valid effective field theory all the way up the the Planck 
scale, it cannot be complete as it contradicts to a number of observations. We would like 
to use this opportunity to underline once more that the confirmed observational signals 
in favor of physics beyond the Standard Model which were not discussed in this paper 
(neutrino masses and oscillations, dark matter and baryon asymmetry of the Universe) can 
be associated with new physics below the electroweak scale, for reviews see [60, 61] and 
references therein. 7 The minimal model — uMSM, contains, in addition to the SM particles, 
three relatively light singlet Majorana fermions. These fermions could be responsible for 
neutrino masses, dark matter and baryon asymmetry of the Universe. The z/MSM predicts 
that the LHC will continue to confirm the Standard Model and see no deviations from it. 
At the same time, new experiments at the high-intensity frontier, discussed in [64], may be 
needed to uncover the new physics below the Fermi scale. In addition, new observations 
in astrophysics, discussed in [61], may shed light to the nature of Dark Matter. As the 
running of couplings in the i/MSM coincides with that in the SM, all results of the present 
paper are equally applicable to the z/MSM. 
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4.1 Note added 

After our paper was submitted to the electronic preprint archive arXiv (on May 13) a 
number of events happened, which require its update. First, extra corrections to the 
matching procedure at low energy scale, not computed in our work, were found in [112] 
(May 29). 8 Ref. [112] finds agreement with our results on the aa s order in the Higgs self- 
coupling constant. In addition, this paper computed a part of the 0(a 2 ) corrections to the 
top Yukawa and Higgs coupling constants in the "gauge-less" limit of the Standard model, 
i.e. the two-loop terms containing the top Yukawa and scalar self-coupling were accounted 
for. The overall effect of these terms happened to be quite small. The corrections shift 
the benchmark Higgs mass up by 0.2 GeV, and reduce the sensitivity of the results to the 
normalisation point from 1.2 GeV to 0.8 GeV, decreasing somewhat the theoretical error- 
bars, see Fig. 3. 

Second, the discovery of the Higgs-like resonance was announced at CERN by ATLAS 
and CMS collaborations [66, 67]. According to CMS, 

M H = 125.3 ± 0.4(stat) ± 0.5(syst) GeV, (4.1) 

while ATLAS gives a slightly higher value 

M H = 126 ± 0.4(stat) ± 0.4(syst) GeV. (4.2) 

And, finally, an updated results on the mass of the top quark were announced at 
ICHEP 2012 (July 9), see [38, 68]. The combination of the Tevatron results reads: 

m t = 173.2 ± 0.6(stat) ± 0.8(syst) GeV = 173.2 ± 1.0 GeV, (4.3) 

whereas the present LHC value is 

m t = 173.3 ± 0.5(stat) ± 1.3(syst) GeV = 173.3 ± 1.4 GeV. (4.4) 

The central values of the top mass are somewhat higher (by 0.3-0.4 GeV) than those which 
were given by the Particle Data Group at the time we were writing our paper. 

In Fig. 9, which is an update of Fig. 4, we show the changes due to the experimental 
shift of the top mass (we take the Tevatron result as having smaller errors) and due to the 
additional two-loop corrections found in [112] for the ATLAS and CMS values of the Higgs 
mass. The value of the benchmark Higgs mass is inside 2a contours on the mt,a s plane 
for the Atlas and CMS values of Mjj, entering the la contours if theoretical uncertainties 

Preliminary results of Ref. [112] were presented in Ref. [111]. However, the 0(aa s ) threshold correction 
to the Higgs self-coupling is not specified in Ref. [Ill], but is presented here for the first time. Moreover, 
the conclusions reached in Ref. [Ill] disagree with ours. 
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Figure 9. Update of Fig. 4. The values of the strong coupling constant a s and top mass M t 
corresponding to a minimal Higgs boson mass M m ' ln coinciding with the recent CMS and ATLAS 
values 125.3 GeV and 126 GeV, and to the Higgs mass values shifted by approximately la ~ 0.6 GeV. 
The latest Tevatron 68% and 95% experimentally allowed regions for a s and M t are given by shaded 
areas. The dashed (dotted) lines correspond to 1 GeV (2 GeV) uncertainty in the M m i n theoretical 
determination. The red lines in the center correspond to the expected precision from a e + e~ collider. 



are incorporated. The red curves show the shrinking the size of uncertainties in mt and a s 
(1 and 2a) which can be achieved with an e + e~ collider operated as a in tt factory on a 
lepton collider (this estimate is taken from [69] and [70] discussing the ILC physics). 

It is important also to note, that it is difficult to determine which renormalization 
scheme corresponds to the numerical value of mt quoted by experiments. The mt determi- 
nations by the Tevatron and LHC collaborations [68, 71] are based on Monte Carlo event 
generators implemented with LO hard-scattering matrix elements. Although it is plausible 
and likely that the experimental mt values thus extracted are close to the pole mass of 
the top quark, this is by no means guaranteed with the due theoretical rigor [72]. In fact, 
rigorous determinations of the top-quark pole mass from total production cross sections 
yield somewhat smaller values, albeit with larger errors [30, 31, 113]. For the time being, 
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we adopt the working hypothesis that the experimental values of rut [68, 71] correspond to 
the pole mass, bearing in mind that this is probably subject to change once a proper NLO 
treatment of the resonating top-quark propagators is implemented. 

Thus, all our considerations remain in force and call for further improvement of the 
theoretical computations, which should account for all 0(a 2 ) corrections to the mapping 
procedure, and for a construction of a ti factory to pin down a s and rnt- A recent paper 
[113], where uncertainties in determination mt and a s have been analyzed, reached exactly 
the same conclusion. 



A 0(aa s ) electroweak corrections to the top Yukawa and Higgs self cou- 
plings in Standard Model 

The evaluation of radiative corrections to the relations between MS parameters (coupling 
constants) and masses of particles includes two steps: evaluation of radiative corrections 
between the Fermi constant Gp and its MS counterpart [73] (see [74-76] for recent reviews) 
and the evaluation of the radiative corrections between MS and pole masses. 

The one-loop electroweak corrections O(a) to the relation between the self-coupling 
X(fi 2 ) and the pole mass of the Higgs boson was obtained in [77] and to the relation between 
the Yukawa coupling yt and the pole mass of top quark was found in [32] . The corresponding 
ingredients for the 2-loop mixed electroweak-QCD corrections were evaluated in [65, 78- 
82], but has never been assembled. We performed independent (re) calculations of all 0(a) 
and 0(aa s ) contributions. In the following we will denote the on-shell masses by capital 
M and the MS masses by lowercase m. 



A.l 0(aa s ) corrections to the relation between on-shell and MS Fermi con- 
stant 

The relation between the Fermi coupling constant and the bare parameters is as follows 
[73]: 

where ARq includes unrenormalized electroweak corrections and go,m 2 Vo are the SU(2) 
coupling constant and the bare W boson mass (see for details [74-76]). After performing 
MS renormalization this relation has the following form: 

qp = + AGpa + AGpaas + ...). ( A.2) 

where on the r.h.s. all masses and coupling constants are taken in the MS renormalization 
scheme. The one-loop coefficient, Ag* f . q , is known from [73] and for N c = 3, Cp = 4/3 
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and m,b = has the following form: 
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Here, sin 2 9\y is defined in the MS scheme as 



sin^ 6 W = sin 2 9$ s (n 



1 



(A.4) 



g 2 {^ 2 ) + g' 2 (n 2 ) ' m|(/i 2 ) ' 

where g'([J> 2 ) and g(fi 2 ) are the U(l) and SU(2) MS gauge coupling constants, respectively. 
The matching conditions between the MS parameter, defined by Eq. (A.4), and its on-shell 
version, [73], follows from identification 
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where Mz and M\y are the pole masses of the gauge bosons (see detailed discussion in 
[83-85]). The evaluation of the mixed QCD-EW coefficient, Ag FjQ , Qs , is reduced to the 
evaluation of the 0(aa s ) corrections to the W boson self-energy at zero momenta transfer 
and may be written in the following way [86-90]: 
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and Z\Y^aa s we used the results 9 of [85] and uf is the number of 
fermion families (tif is equal to 3 in the SM). 
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Using the fact, that Gf is RG invariant, i.e. ^i 2 -^Gf = 0, the /x-dependent terms in 
Eq. (A. 6) can be evaluated explicitly from the one-loop correction and explicit knowledge of 
anomalous dimension jg f - As was shown in [83-85, 91, 92], the anomalous dimension jg f 
can be extracted (i) via the beta-function (3\ of the scalar self-coupling and the anomalous 
dimension of the mass parameter m 2 (in unbroken phase) or (ii) via the /3-function of the 
SU(2) gauge coupling g and the anomalous dimension of the W boson (in broken phase): 



1G F = j" 

Eq. (A. 2) can be written as 
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Collecting all terms in Eq. (A. 10) we get 
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At the end of this section we again point out that the anomalous dimension of the 
vacuum expectation value v 2 (fj, 2 ) = 1/(\^2Gf{h 2 )) within the diagram technique is defined 
by Eq. (A.7) and it is not equal to the anomalous dimension of the scalar field as in 
the effective potential approach [93]. Another important property of Eq. (A.7) is the 
appearance of an inverse power of the coupling constant A due to the explicit inclusion of 
the tadpole contribution. As consequence, the limit of zero Higgs mass, m 2 H = 0, does not 
exist within the perturbative approach. The importance of the inclusion of the tadpole 
contribution to restore gauge invariance of on-shell counterterms was recognized a long 
time ago [94] and was explicitly included in the one-loop electroweak corrections to the 
matching conditions [32, 77]. The RG equations for the mass parameters were discussed 
in [83-85, 91, 92]. 
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A. 2 0(aa s ) corrections to the relation between the MS and pole masses of the 
top quark 

The detailed discussion and explicit evaluation 10 have been presented in [82] (the results of 
[82] were also used for analysis of convergence of series representation of the set of Feynman 
Diagrams evaluated in [96, 97]). For our analysis is enough to write the following symbolic 
relation between the MS and pole masses of the top quark: 
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where a a and a aoia are defined by Eq.(5.54) or Eq.(5.57) of [82]. 

The pure QCD corrections can be found in [24-26] (only the value of a aa (M t ) is given 
there, but the expression for other \x values can be readily reconstructed from the beta 
functions). 

A. 3 0(aa s ) corrections to the relation between the MS and pole masses of the 
Higgs boson 

At the two- loop level the relation between the pole and MS masses is defined as follows: 
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where the sum runs over all species of particles, gj = a, g s , ((Jg^o) and (£m| ) are 
the one-loop counterterms for the charges and physical masses in the MS scheme and after 
differentiation we put all parameters equal to their on-shell values. The derivatives in 
Eq. (A.14) correspond to the subtraction of sub-divergences. The genuine two-loop mass 

masses of the Higgs boson has the form 
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where g is the SU(2) MS renormalized coupling constant. 



There is typo in Eq. (4.46) of [82]: the common factor Cj was lost. The correct result is 
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However, all plots, the Eq. (5.57) and the Maple program [95] are correct. 
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The exact analytical result for the 0(aa s ) two- loop quark contribution to the Higgs- 
boson self-energy was calculated in [81, 88]. The bare two loop mixed QCD-EW contribu- 
tion (with explicit inclusion of the tadpole) for the quark with mass m q reads: 
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In accordance with Eq. (A. 14), the coefficient A m 2 )Q!aa) g of order 0(aa s ) relating the pole 
and MS masses of the Higgs boson, s p — m 2 H , can be written as 
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As result of our calculation we find: 
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The contributions of other quarks with non-zero mass are additive. Exploring the e expan- 
sion for the master integral Jo qq from [99], we have for t-quark contribution (q = t): 
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where 
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and we have introduced the two functions F(y) and G(y) (see also [80, 81]) defined as 11 



F(y) = 3 [Li 3 (y) + 2Li 3 (-y)] - 21ny [Li 2 (y)+2Li 2 (-y)] 

-*ln 2 y[ln(l-y) + 21n(l+y)], 
G(y) = [Li 2 (y)+2Li 2 (-y)]+lny[ln(l-y) + 21n(l+y)] 



and 
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In Eq. (A.21) we explicitly factorized the RG logarithms, C H ' aa and C^' aa , which may be 
calculated also from the one-loop result and the mass anomalous dimensions (see [91, 92] 
for the general case). From the parametrization 
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n We cross checked, that Eq. (A.21) minus tadpole contribution coincides with results of Ref. [80, 88-90] 
after the following substitutions: 
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11 

where L = In and using the fact that pole mass is RG invariant, we have: 
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In terms of the variable y, defined by Eq. (A. 22), the final result reads: 
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A. 4 0(aa s ) corrections to the top Yukawa and Higgs self couplings 



The relation between the top Yukawa (Higgs) coupling and the Fermi constant G F is 
obtained from Eqs. (A. 2), (A. 13) and (A. 26) as: 
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and the relation between the Higgs coupling constant A = fygwin used in [77] and the 
parametrization of [82-85] follows from the comparison of the RG functions: /isirlin = 

Ajegerlehner(M 

The 0(aa s ) result for the top- Yukawa coupling reads (see Eq. (21) in [65] and [100]) 
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where ax are defined in Eq. (A. 13). The 0(aa s ) result for the Higgs coupling is 
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and the sum runs over all quarks. 

For completeness, we present also the explicit expressions for the derivatives (for iV c 
3, Cp = 4/3 and m& = 0): 
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Terms of the order 0(a), 0(a s ) in Eq. (A. 34) and Eq. (A.35) correspond to [32] and 
[77], respectively. Terms of the order 0(aa s ) in Eq. (A. 34) and Eq. (A.35) are the 
mixed electroweak-QCD corrections and Ag Pi(Ms , cr aois , A m 2 aoig are defined by Eq. (A. 6), 
Eq. (A.21), and Eq.(5.54) or Eq.(5.57) of [82]. 

For completeness we present also the the coefficient A m 2^ a . We divide all corrections 
into bosonic (diagrams without any fermions) and fermionic (diagrams exhibiting a fermion 
loop) ones: A m ^ a = ^m 2 H (A m ^ Qjboson + A m2 ^ fermion ) . Using the notations of 
[83-85] we may write the one-loop corrections in the following form 
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where (see Eq. (E.6) in [99]) 
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All results are collected in the Maple code of Ref . [95] . 
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B Beta functions 

Two loop SM beta functions above the top mass are collected in [10] (see [101-110] for 
original works). The three loop beta functions can be read off [19, 20]. 

Below the top mass the one loop beta functions for the gauge couplings were used to 
evolve the PDG values from Mz to Mf. For example, for the a(/i) 
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The higher loop corrections are not important numerically for the electroweak constants 
for the small energy range between Mz and Mf. 

For the strong coupling a s = g|/(4-7r) the RG equation up to order 0(a s 3 ) is used 
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and Nf = 5 is the number of flavors below the top quark. Strictly speaking, the value of 
ots(Mt) obtained from this equation should be also shifted to the 6-quark value by 



11 



as,N f =6{M t ) = a StNf= s(Mt) - ^^a SNf=5 {M t ), 



(B.3) 



but this introduces a negligible effect (< 0.1 GeV) for the Higgs mass. 

In all the formulas of the Appendix A we use the values of a and a s at the matching 
scale \x. 
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Figure 10. Contributions to the top Yukawa constant from QCD up to 2 loops, up to 3 loops, 
QCD and 1 loop EW corrections O(a), and QCD with 0(a) +0(aa s ) corrections. One parameter 
is vrying, the two others are chosen from M t = 172.9 GeV, a s — 0.1184, Mh — 125 GeV. 
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